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Abstract 



^ . We analyze the Farey spin chain, a one dimensional spin system with 

O I effective interaction decaying like the squared inverse distance. Using 

^ ■ a polymer model technique, we show that when the temperature is de- 

creased below the (single) critical temperature Tc = ^ , the magnetization 
jumps from zero to one. 



1 Introduction 



Can a magnet keep its full mean magnetization (m) = 1 up to the Curie tem- 
perature Tc and then loose it at one stroke? Definitely such a property would 
be different from the usual situation, where (m) continuously decreases to zero 
(though not being differentiable at T^), or jumps discontinuously by an amount 
^ ■ strictly less than the saturation value. 

■ It has been proven [0, ^ ^, ^, |T3|, |T5| that certain spin chains of long range 



ferromagnetic interaction exhibit a discontinuity of (m) at Tc, jumping from a 
value in the interval (0, 1) to zero. 

In one dimension such a phenomenon can only occur if the effective interac- 
tion between spins of distance d decays at most like d^"^, since there cannot be 
a phase transition for a decay rate of d~" if a > 2. 
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However these examples do not exactly provide a positive answer to the 
question posed initially, since the jump of (m) at Tc is strictly smaller than 
one. Indeed for non-zero temperatures a mean magnetization (m) = ±1 is 
only possible if the forces between the spins become so strong that one may 
legitimately ask whether this could endanger the existence of a thermodynamic 
limit. 

In the present paper we show the contrary by considering the example of the 
Farey fraction spin chain. Similar to the abstract polymer model formalism, 
is introduced in Sect. 2. Since the limit free energy coincides with the one 
of the number-theoretical spin chain [Theorem^), the single phase transition 
(nonanalyticity of the free energy density) is situated at inverse temperature 
[3 = 2. 

In Sections 5 resp. 6 we consider the mean square magnetization in the 
regimes below resp. above the temperature. Whereas (m?) (/?) = 1 for low 
temperatures {Theorem 0), this quantity vanishes above Tc {Theorem 0). 

We conjecture, and plan to prove, that the spin chain has exactly two extremal 
Gibbs measures in its low temperature phase, and one above Tc . 

We also invoke a polymer model technique similar to the one developed in 
[H to estimate the strength of the interaction. 



2 The Model 



in IjT^ the so-called number-theoretical spin chain was introduced, whose low- 
temperature partition function equals a quotient of Riemann zeta functions. In a 
series of subsequent papers (see for a survey) this model was then analyzed 
further. In particular it was shown in that a phase transition with a jump of 
m from one to zero occurs at = |. 

The number-theoretical spin chain shows an asymptotic decay of interactions 
which is exactly of the form d^"^, and the limit free energy density exists. The 
main motivation of its study lies in its connection with number theory, and more 
specifically in the hope that its ferromagnetic character together with a version 
of the Lee- Yang theorem could shed a light on the location of the zeroes of the 
Riemann zeta function. 

^From the statistical mechanics point of view it should, however, be said 
that it lacks the strict symmetries usually encountered in ferromagnets. It is 
neither fully translation invariant nor invariant under spin reversal, although both 



2 



symmetries are asymptotically present in the bulk. 

In the so-called Farey fraction spin chain was introduced as a spin system 
of statistical mechanics related to the Farey fractions in number theory. As we 
shall state below, this chain, which has strong relations with the one mentioned 
above, but a less direct number-theoretical interpretation of its partition function, 
has all relevant symmetries. 

The definition of the Farey chain in was based on functions 

Mfc : Gfe SL(2,Z) (A;>0) 

on the additive group Gk '■= {0, Ijii' - '''}, inductively defined by setting Mq : = 
(J 5) and for A; > 1 

Mfc(a) := A'-'^^B'^'-Mk-iiai, . . . , cxk-i) (a G G^), (1) 
with A := (l °) and B := = (J }). The function 

Ek := In(Tfc) with Tk := Trace(Mfc) : Gfc ^ N 
was interpreted as the energy function of a spin chain with k spins with values 

CTi, . . . , (Tfc. 

Then by discrete Fourier transformation 

(J^Mt) := 2-' J2 /(^) • i-^r" (t e G,) (2) 

the energy function has the form 

Ek{(r) = -J2Mt){-ir' (fxGGfc) 

with the so-called interaction coefficients 

Jk{t) := -{TkEk){t) (tGGfe). 

The 'lattice gas' spin values cTj = 0, 1 are used here for convenience. The mean 
magnetization 



irtk : = 




however, is defined using the spin values Si{a) := {—lY' G {±1}- 
The Farey spin chain has the following symmetries: 
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1. When one interprets {1, .... A;} as a system of representatives of the residue 
class ring Z/ZcZ = {/ + /cZ | / G Z}, then by cyclicity of the trace the 
energy function is invariant under the shift 

Sk : Gk ^ Gk , Sk{(T)i :— ai-i (3) 

on the configuration space of the chain. So the interaction is translation- 
invariant, too ( JkO Sk — Jk)- 

2. Since AP = PB for P := Q J) 

Mk{ak, . . . , (Ti) = PMfe(ai, . . . , akfP. 
This implies the mirror symmetry 

Ek{ak, . . . , (Ti) = Ek{ai, . . . ,ak) 

and a similar relation for the interaction coefficients. 

3. Finally we notice that by 2) and transposition invariance of the trace 

Ek{l - a) ^ Ek{a) for 1 - := (1 - (7i, . . . , 1 - (7^) 
so that 

k 

Jfe(i) = for ^ijodd. 

i=l 

By 3) we need only consider t in the even subgroup 

Gl-= {t^GklEliti even}. 

Jfe(O) < 0, since this is the negative mean of the (positive) energy function 
Ek- Note that this is the only interaction coefficient which does not influence 
the Gibbs measure. 
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3 A Polymer Model Interpretation 



The notion of polymer models grew out of an abstraction of situations like the 
one encountered in the low temperature expansion of the Ising model. There one 
may decompose contours X into non-intersecting cycles "fi { X = (71, . . . ,7/)), 
and express the Boltzmann factor of the spin configuration in terms of products 
of activities attributed to these cycles (the activity zi^ji) of a cycle equals the 
exponential of its length, multiplied with minus the inverse temperature). 

In an abstract setting (see, e.g., Gallavotti, Martin-Lof and Miracle-Sole 0, 
Glimm and Jaffe [[7|] and Simon |]T^) one starts with a set P (which we assume 
here to be finite), whose elements are called polymers. Two given polymers 
71,72 G P may or may not overlap (be incompatible). Incompatibility is assumed 
to be a reflexive and symmetric relation on P. 

Thus one may associate to a I- polymer X := (71, ... ,7;) G an undirected 
graph G{X) = {V{X),E{X)) with vertex set V{X) := {1,...,/}, vertices 

j being connected by the edge {i,j} E E{X) if 7^ and 7^ are incompatible. 
Accordingly the /-polymer X is called connected if G{X) is path-connected and 
disconnected if it has no edges ( E{X) = 0). 

The corresponding subsets of are called resp. D', with := P° : = 
{0} consisting of a single element. Moreover P°° := [J'^qP^ with the subsets 
■= UZo 3nd := C^- We write |X| := / if X G PK 

Statistical weights or activities 2; : P — > C of the polymers are multiplied to 
give the activities := Y[\=i ^ili) °f /-polymers X . A system of statistical 
mechanics is called polymer model if its partition function Z has the form 



Air W 



Then, up a normalization factor, the free energy is given by 

with n{X) := n+{X) — n_(X), n±{X) being the number of subgraphs of G{X) 
connecting the vertices of G{X) with an even resp. odd number of edges (see 
Gallavotti et al. |||]). It is known that (see, e.g. Prop. 20.3.5 of 10]) that 

(-l)l^l-^n(G) > 0. (6) 
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This also follows from the deletion-contraction property 



n{G) 



n{G') - n{G") 



(7) 



where G' is obtained from G by deleting an edge and G" is the graph which 
arises by contracting the same edge of G (see, e.g. Read ||12|| ). 
In the present context of a chain with k spins we use 

• the set Pk of 1 + k{k — 1) polymers given by 



• We map the polymers 7 G to group elements 7 G by setting p := 
and pi^r := 5i + 5^, where for i G Z/ZcZ the group element 5i has the form 
5i{l) = 1 if / = 2 and zero otherwise. 

This map induces a map 



from the set of multi-polymers to G^. 

• The support of our polymers is given by supp(p) := 'L/k'L and 

supp(pi,r) := {/, / + 1, . . . , r — 1, r} C Z/A;Z. 

Note that the polymer p^,; is different from pi^r, although the group el- 
ements pi^r and pr,i coincide, and for /c = 2 the supports supp(p) = 
supp(pi,2) = supp(p2,i)- 

The polymers 7 and 7' are called overlapping or incompatible if 
supp(7) n supp(7') 7^ 0. 

• We attribute to the polymers the activities 



Pk := {p} U {pi,. \ l^re Z/kZ}. 



X = (71,..., 7i) ^X:=^7i 



i=l 



Z 



(p) ■= 3-N"PP(p)l = s-k and z{pi^r) ■= -S'I^^pp^?''''-)!. 



(8) 
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Every group element t G allows for exactly two representations t = X hy 
disjoint multi-polymers X E D^. 

Lemma 1 The Fourier transform jk := J-'kTk can be written as 

Mt) = {-^' J2 <x) (teG,). (9) 



XeD^:X=t 



Proof. 

• For t odd both sides are zero. 

• For t = we perform the sum to obtain 



j^(0) = 2-^ Trace(Mfc(cT)) = 2-^Trace(^^) 



o-eGfc 

with 5* := A + B = [H), which has eigenvalues one and three. So 
j;,(0) = {3'' + l)/2'=. On the other hand, the r.h.s of (D is of the form 

(|)'=(^(0) + ^(p)) = (f)^l + 3-^). 

For t E Gl \ {0} we assume without loss of generality, using cyclicity of 
the trace, that 

with 0„ = (0, . . . , 0) e Gm- Then with D:= A- B= {\ "J) 
jfc(t) = 2-'^Trace . . . DS'^^'^D) . 

Now 

DS^'D = S"" - {3"^ + 1)1 (10) 
commutes with S, and D'^ = —1 so that 



with Ami := X]r=i "^2«-i and Am2 := X^ILi'^si- Thus using ( pD]) we 
arrive at 

This equals the r.h.s. of (H). □ 

4 Comparison with the Number-Theoretical Spin 
Chain 

The number-theoretical spin chain of length k has the canonical energy function 

:= In(h^) with : Gfc ^ N 
inductively defined by 

h^:=l, h^+i(a,afc+i) :=h^(a) + afe+ih^(l-a), (a G Gfe).(ll) 
It turns out to be useful to consider the grand canonical energy functions 

H^iGfc^M with H^(a) :=H^+,(a,l) (a G G^), 
too, which is the logarithm of 

h^Gfc^N , h^(a):=h^+i(a,l) = h^(a) + h^(l-a). (12) 

Namely in polymer model techniques were applied to estimate the grand 
canonical interaction, := — JF^H^. These were applied to the subset 

Pk := {pi,r ePk\l<r} (13) 

of polymers (where the inequality < in "L/kZ is understood as the one for the 
representatives in {1, . . . , k}). This is the set of polymers which contribute to 
the thermodynamic limit. 
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For t G Gfc \ {0} the resulting formula 



J^it) = -5t,o ■ (ln(2) + Hn(3/2)) - ^^"^ 



X=t 



\x\\ 



(14) 



for these grand canonical interaction coefficients contains only nonnegative terms. 
This follows from (|^) and the fact that all activities (|^) of polymers in i\ are 
negative. Similarly the canonical interaction of the number-theoretical spin chain 
was shown to be ferromagnetic. 

Lemma 2 The Farey interaction coefficients Jk{t) = —^kEk{t) can be written 
as 



Jk{t) = ~St,okH3/2) - ^^"^ 



|X|! 



x=t 



Proof. Since jk ■= J^kTk, we have Ek = 2^J^kjk and 



-2-- 1„ 



u-t 



-5t,o-A;ln(3/2)-2-'= J] In 

<TGGfe 



zX 
IXI! 

xgd?° ' ' 



-1) 



a-t 



where the redefined single-polymer activities z^ji^^), E Pk are given by 

z.(7) := z^{^) ■ {-ir\ 
that is z„{p) = z{p) and z„{pi^r) = z{pi^r) ■ {-ly'^"'' ■ By (|) we get 
Jk{t)+6t,o-k\n{3/2) 



cr 

■ E 



(xGGfc X&C^ 

n{X) . 

A/ 



\X\\ 



\x\\ - 



crGGi 



(15) 



(16) 



(17) 
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using the identity E<xeG,(-l)"" = '^^s,o- □ 

Although formula (|T^) looks very similar to (p^l), the sum is over all connected 
multipolymers based on the full set Pk of polymers, instead of the subset (|13|). 
Therefore not all terms in that sum are positive. By (|]) and (|) the negative 
terms are precisely the ones containing an odd number of copies of the polymer 
p. Thus the positivity of the interaction for finite k does not follow immediately. 



5 The Free Energy 

Theorem 3 Tlie limit free energy density 

F{13) := lim Fu{(3) of Fu{(3) := ^ ln(Zfc(/3)) {(3 > 0) 

witfi k-spin partition function Zk{P) ■= JZaGG^ 6xp(— exists and equals 
the one of the number-theoretical spin chain. 

Proof. We use the canonical and grand canonical ensembles as bounds for F^. 
1) Since the entries of the matrices Mfc(cr) are non-negative, an upper bound 

for Tfc = Trace(Mfc) is given by Trace ((J J)Mfc). But for a e Gk 

Trace ((; S)Mfc(a)) =h^(a) and Trace ((J ;)Mfc(a)) = h^(l - a), (18) 

since both sides equal one for k = 0, and 

Trace ((; l)Mk{a)) 
= Trace ((} ^A'-^" B^" Mk-l{a^, ak-i)) 
= Trace (((J J) + a^J J)) Mfc_i(ai, . . . , 

= h^_i(cri, ■ ■ ■ , o-fe-i) + o-fch^-i(l - cxi, . . . , 1 - CTfe.i) = h^(o-) 



and similar for the second identity in ( ^81) . Adding these identities and, using 
Def. dll), shows that 



Tk < (ke H 







2) To derive a lower bound for Tk, we notice that for a E Gk-i 
Trace ((J J)M,(0,a)) = h^.ia) and Trace ((J ;)Mfc(0,a)) = hti(l - a) 
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since both sides equal one for k = 1, and for a G Gfc_i 

Trace ((; l)Mk{0,a)) 
= Trace ((J J) A4_i(0, ai, . . . , 

= Trace (((;;)+ cx^-i Q J)) M,_i(0, ai, . . . , a,_2)) 
= h^_2(cri, • • • , crfe_2) + o-fc_ih^_2(l - c^i, . . . , 1 - o-fc_2) = h^-i(o-) 
and similar for the second identity. Thus 

Tfc(0, a) = Tfe(l, 1 - a) > h^iM (A; G N). (19) 
3) Since the (grand) canonical free energies are given by 



In ( E.GG, ( 



these two inequalities imply 

The canonical and grand canonical ensembles have the same limit free energy, 
since 

F^<F^<F,^+^^^^t±ll. (20) 

So the limit free energy F of the Farey chain coincides with the one of the 

number-theoretical spin chain 

The lower inequality in (|20| ) follows from (0), the upper inequality from 
< + 2) ■ h^, which is a consequence of (0) and the relation 

h^(l-^)<(A: + l)-h^(a) (aeG,) 

(which follows from Def. (|ll]) by induction). □ 
The same conclusion was reached in by a different method. 

Corollary 4 The Farey spin chain has exactly one phase transition, at (3 = 2. 

Proof. This follows from the corresponding statement in for the number- 
theoretical spin chain. □ 
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6 Low Temperature Magnetization 



Due to the invariance of the energy function Ej. w.r.t. spin flips, the mean 
magnetization has expectation zero. However, the long-distance correlations 
are measured by the square of that variable. 

Theorem 5 In the low temperature phase (3 > 2 

lim(m^>J/3) = l. 

Proof. We complement estimate (|^) by 



which follows inductively from Def. (|i|) by noticing that for such a both off- 
diagonal entries are > 1. We thus have 



Tk{cj) > k 



(aeGfc,(0,...,0)^a^(l,... 



1)), 



oo 




(21) 



n=l 



with 



(^k{n) := \{a e Gk \ max(h^(cr), k + 1) 
It is known [jTO|] that 



n}\. 



can be written in the form 



oo 



n=l 



with (fk{n) < fin) and (pk{n) = <f{n) for n < A; + 1, 



ip{n) := \{i e {!,... ,n} \ gcd(i, n) = 1}\ 



being Euler's (/^-function. Thus for Re(/3) > 1 
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Substituting the identity 

oo k 
n=l n=l 

for the r.h.s. of (^), we thus get 

k 

Urn \Zk{f3) -2 -2-^1 < lim |(A; + ■ y (^(n) 

fc^oo fc— >oo ' ^ 

n=l 

fc 

< lim + ■ Vn = 0, 

fc— +CO ' 

n=l 

SO that Mmk^oo Zk{f3) = 2-2'^ for Re{(3) > 2. Using < < 1 and 
m2.((0, . . . , 0)) = m|((l, . . . , 1)) = 1, we thus get limfc^o, (mf,), (/?) = 1. □ 

This extends the same conclusion, reached by a different argument for /? > 3 f^. 

For the ferromagnetic spin chain the limit mean magnetization (m) := 
liiiifc^oo ("^fc)fc equals 1 for the canonical ensemble and (3 > 2, whereas it van- 
ishes in the high temperature region [If. 

For the grand canonical ensemble, as for the Farey ensemble, (m) vanishes 
identically, since the interaction is even. Of course this does not say much about 
the actual structure of the extremal Gibbs states. 



7 High Temperature Demagnetization 

Now we consider the mean magnetization in the high temperature regime (3 <2. 
To show that the expectation of the square vanishes, we need a correlation 
inequality. So consider for n G N the configuration r G G„+2 with spins ri : = 

Tn+2 := , := 1 for 2 < / < n + 1 , and the event 

81 := {a G Gk+n+2 I ff; = for 1 < / < n + 2} 

of an initial string of n adjacent 1-spins enclosed by 0-spins. 

Due to the long range character of the interaction, one might think that, 
given 8]^, the ferromagnetic interaction would tend to align the other spins in 
the 1-direction (equal to r2 = . . . = r„+i), at least if n is large. 
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This would mean a negative conditional expectation of Sj = (—1)^' for i e 
{n + 3, . . . ,n + k + 2}. Because of the dominance of the multi-body interactions 
this is, however, not the case. 

In fact, the non-inverted spins (—1)"^^ = (—1)^"+^ = 1 tend to produce an 
anti-ferromagnetic effective coupling between the spins in the regions 2, . . . ,n+l 
and n + 3, . . . ,n + k + 2 they separate: 

Proposition 6 For A C {n + 3, . . . ,n + k + 2} and /S > 

(sA|^nfc+„+2(/3)>0 with SA=l[s,, (22) 

(/ \£)i denoting the expectation of f : G; — R , conditioned by the event £ . 
Proof. We set 

T,":Gfc^N , r,"(a):=Tfc+„+2(r,a) 

and := \n{T^) . We first prove 

(^fci?fc)(t)<0 (tGGA{0}), (23) 

using a polymer technique similar to the one above. Namely we redefine the set 
Pk of polymers by 



l<l<r<k ' 



and map them to the group elements := p^ := Sm G G^ resp. pi^r '■ = 
5i + 5r E Gfc. Depending upon the length n of the 1-substring, the polymer 
activities are given by 



o-|supp(p^)| o-|supp{p«)| 



'Lr) 



where 

supp(p^) := {1, . . . , m}, supp(p^) := {m, . . . , k} and supp(pi,^) :={/,..., r}. 

Polymers with intersecting supports and the polymers p^ , p^ are mutually 
incompatible. We now claim that in analogy with (^ the Fourier transform 
:= J^kTJ} can be written as 

m=2in + l){lY (^^G,). (24) 

XeD?^:X=t 
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To show this, we write t = (ti, ... ,4) uniquely in the form 



j^(t) = 2-'=Trace {NS"''DS'''W . . . S""-) 
with N := AB^'A = {n + 1) ■ (H) + D . 
• If M is odd, then 

j^(t) = {n + l)- (-l)("-=')/22-*^Trace ((} S^'^WS'^'^W) 

with Ami := Z]i=i ^^^^ "^2i-i and Amg := Zli=i^^^^ "^2i- So 
= 2{n + 1) ■ (_i)(«-i)/22-fe3Ami 



with Ami := Zirii "^2i-i and Ams := I]rii"^2i- So using (|0D 




(mi > 0) 



so that 



2(^+1)- (I)' n (-3-'"^'~')- 



• If M is even, then 
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In both cases this coincides with the r.h.s. of Similar as in Lemma 0, we 

get 



x=t 



from which ( ^3] ) follows, using and the negativity of all polymer activities. 
Now 



SO that ( ^21) is a consequence of the first GKS inequality for ferromagnets. □ 
Theorem 7 In the high temperature phase < /5 < 2 

lim (m^) {(3) = 0. 

Proof. Since by translation invariance {rrig)^ = ^ Yl^j=i {■^i-^j)g' suffices to 
show that for £ > there is a uniform correlation estimate of the form 

I {sis,)^ m <e {je iUe), . . .,g-j,(e)}). (25) 
We consider the family < £^ol'„_2 L=i,. .nmax of events 

I ^ ) 1=1,. ..,n 

^;-n-2 := '5,-'^^;-„-2 C G; with C^^ := {a G GJ a, = l}, 
using the shift map on Gg. As these events are disjoint, 

E PA.(^;i-2) = n,, ( u ) < P/3,,(Gj) = i (26) 

for the Gibbsian probability P/3,fe(^) := E^e^ ^"^^"^"V^fcl/?) of an event ^. On 
the other hand if /3 < 2, for e > there is a n^g^Je) with 



EEP/^-^(^;-n-2)>i(l-^) (27) 



n=l 1=1 
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for all large g. This property, which is specific to the high-temperature region, 
can be proved as follows. We note that the thermodynamic limit of the internal 
energy 

U:=\un^Uk with f/^ := (i^fc)^ 

exists and equals U{(3) = By concavity and analyticity of (3 ^ PF{f3), 

and by F{f3) = for /3 > 2 we conclude that 

UiP) > iP< 2). 

This implies a positive limit density of spin flips between neighbouring spins and 
thus the existence of an n^g^^^e) meeting (P7|). 
Since by spin inversion symmetry 

for Cg := Gl \ U„,;^;i„_2, by (|6D, (|7D and Proposition | 

{sis,)g < 2Pft,(C,) < e. 
Together with a converse estimate this proves (|25|) . □ 
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